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Abstract 



We present a generalized circle theorem which includes the Lee- Yang the- 
orem for symmetric transitions as a special case. It is found that zeros of the 
partition function can be written in terms of discontinuities in the derivatives 
of the free energy. For asymmetric transitions, the locus of the zeros is tan- 
gent to the unit circle at the positive real axis in the thermodynamic limit. 
For finite-size systems, they lie off the unit circle if the partition functions 
of the two phases are added up with unequal prefactors. This conclusion is 
substantiated by explicit calculation of zeros of the partition function for the 

Blume-Capel model near and at the triple line at low temperatures. 
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More than three decades have passed since Yang and Lee |J first published their cel- 
ebrated papers on the theory of phase transitions and the circle theorem on the zeros of 
the partition function of ferromagnetic Ising models. Although there have been some devel- 
opments in extending the theorem 0, little is known about what happens to the theorem 
in the case of more general first-order transitions such as asymmetric transitions or tem- 
perature driven transitions. In their original proof of the theorem, Lee and Yang relied 
on particular properties of the coefficients of the polynomial. Furthermore they used an 
analogue of a certain electrostatic problem to relate the discontinuity in the spontaneous 
magnetization to the density of zeroes. In this paper we present a simple proof of the circle 
theorem which can be extended to first-order phase transitions of the most general kind. 
The extended theorem can be used to (1) identify the order of the phase transition, since a 
second-order phase transition manifests drastically different characteristics in the complex 
partition function ||; (2) pinpoint the transition point very accurately; and (3) resolve the 
recent controversy over the equal weight versus equal height of the probability distribution 
function at an asymmetric transition [[§-§] among other potential applications. 

The generalized theorem is based on a very simple property that the probability distribu- 
tion for the order parameter M, conjugate to an external ordering field H is doubly peaked 
in the two-phase region. This simple property alone is sufficient to prove the generalized 
circle theorem which states that the locus of zeros of partition function in the complex em- 
plane forms a circle near the positive real axis. Near a triple line where three phases coexist, 
there are three peaks in the distribution function, and this leads to some interesting behav- 
ior of loci of complex zeros. These conclusions are substantiated by explicitly calculating 
zeros of the partition function of the Blume-Capel model near and at the triple line at low 
temperatures. 

Let us consider a partition function Z(H) of a discrete lattice system made of a d- 
dimensional cube of side L as a function of variable H; we seek the complex zeros of Z(H). 
H can be an external magnetic field for a spin system in which case the conjugate order 
parameter would be the magnetization. 
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In the following we will consider only three phase coexistence for the sake of simplicity. 
Generalization to coexistence of more than three phases is straightforward if cumbersome. If 
we label the three phases as A, B and C and the partition function of each phase as Z\(H), 
where A = A, B, C, then we have near a triple point 

z(h)=j:z x (h). (i) 

A 

The fact that the partition function is made of three parts follows from the fact that 
there are two barriers, say at K\ and K 2 in the order parameter M, which separate three 
phases. Therefore we can define Z A (h) = EM=M min e hM Z M , Z B (h) = Em^kI^Zm and 
Z c {h) = EM>K? ax ehMz M, with h = H/K B T = f3H. Z M is a partition function for fixed 
M in case H is the external field and we may take 0>(E), the number of configurations at 
energy E, when H = T and M = E for temperature driven transitions. We can always 
make h — at the coexistence point by redefining Zm appropriately. 

Let M be max(|M min |, |M max |). Introducing the probability distribution function 
p(x) = Z m=x m/Z(H = 0), we can define the moment generating function M.{t) = (e tx ) = 
J2 x e tx p(x) = YliZo( xl )t l /l" Since \\(x l )\\ < 1, M{t) is analytic in the whole complex 
t-plane. Since M(t = hM) = Z(ti)/Z(0), zeros of partition function Z(h) can be ob- 
tained from zeros of Ai(t). By separating the contribution from three phases, we can write 
M(t) = M A (t) + M B {t) + M c (t), where M x (t) = Z x (h = t/M)/Z(0). 

We first consider the case where only two phases A and C coexist. Introducing the 
cumulant generating function defined by ip\(t) = \n(J\4\(t)), we can rewrite Ai(t) = 
2 e ^(*) cosh(-^(t)), where ip = (ip c + V ; a)/2 and $ = (ip c ~ V ; a)/2. Zeros of the M(t) in 
the complex t-plane are simply solutions of following equation, ip(tk) = ±i(l/2 + k)ir, where 
k = 0, 1, 2, 3 • • •. Zeros of partition function in the complex can be calculated from the 
relation h k — t k /M. 

In order to express the zeros of Z(h), in terms of the free energy density and its derivatives 
with respect to h of each phase, we expand ip\ in Taylor series of h as ip\ = ln(A^A(0)) + 
L d EZi h l rf/tt, where ^ = L~ d d l \n(Z x (h))/dh l \ h=0 . Therefore tp(t = hM) = - ln(a)/2 + 
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F(h), where a = M A (0)/M c {0) and F(h) = L d YZ,i hl li/V- wi th li = {if - 7, A )/2. 
Therefore zeros of Z(h) are solutions of the following equations, 

%(F(h)) = ±(1/2 + k)7r = I k (2) 
U(F(h)) = ln(a)/2 = R. (3) 

We can easily obtain zeros Z(h) by inverting the series. Following the tradition, we will 
consider zeros in the z-plane defined by z = e h . The modulus r = \z\ = e^ h > and argument 
9 = ln(z/r)/i = Im{h) of zeros can be found from 

ln(r fc ) = U(h k ) =R + U-R 2 + 4 2 ) + (27 2 2 - 7s) 
x(R 3 - 3RI 2 ) + (-5 72 3 + 5 727 3 - 74) (R A - 6R 2 I 2 k + %) 
+ (14 72 4 - 2l7 2 2 73 + 3 73 2 + 67274 - %)(& - WR 3 P k 
+5Rl£) + (42 72 5 - 84 72 3 73 + 28 72 7 3 2 + 28^74 - 7 73 74 

-77275 + 7 6 )(-^ 6 + l5R 4 P k - IhRHi + It) ■■■ (4) 

e k = %h k ) = i k - 2%M k + (2f 2 - %)(3R 2 i k - il) 

+4(5 72 3 - 5 7 273 + %)(-R 3 Ik + Ml) + (14 72 4 - 2l7 2 2 7s 

+3 73 2 + 67274 - %)(5R% - ioi? 2 4 3 + il) 

+ (-42 72 5 + 84 72 3 7 3 - 28 72 7 3 2 - 287^74 + 7 73 74 + 7 72 7 5 

-7 6 ) (6i? 5 4 - 20,R 3 4 3 + 6i?4 5 ) ■ ■ • (5) 

where R = R/L d fi, I k = I k /L d fi, 7; = iij \xl\ and \i = 71, is one half of the discontinuity 
of the order parameter across the phase boundary 

The angular density of the zeros, g{9) defined by L d g{6) = \/{9 k+ i — 9 k ), can be obtained 
by differentiating the left hand side of (0) with respect to 9 k . We have 

2ng(9) = 2^ - %6 2 + %9 A /12 + 0{L~ d ) (6) 

In the scaling limit, we have r = exp[ln(a)/(2L d /i) + 0(L~ 2d )] and 9 k = ±n(l/2 + k)/(L d n) + 
0(L~ 2d ). This implies that the zeros are distributed uniformly on a circle of varying radius 
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near the positive real axis for finite system if the partition functions of the two phases are 
added up with unequal prefactors. If the prefactors are equal, i.e, a — 1, then zeros form a 
circle of unit radius near the positive real axis in the scaling limit. However if the argument 
, 6k grows large enough zeros deviate from the unit circle unless all 77's with even I vanish. 
However for symmetric transitions we have not only a = 1 but also 7; = for all even I. 
Therefore we have = 1 for all sizes and the angular distribution of the zeros becomes 
exactly Eq. @ without the size dependent correction terms. This is the Lee- Yang theorem. 
We now have the angular density of zeros in a complete power series in 9. 

Let us now consider a case where three phases coexist. Using Fs{h) and 6 defined by 
F B (h) = L d Y,Zi h l lf/l\ and 6 = M B (0)/M c (0), we can modify eqs.(|) and (|) as 

e JKF(*)) (e 8(^)) + ae -»(F(h))) cos(3F(/0) 

= _ he ^{F B {h)) cos^Fs^))), (7) 

e W/0) (e R(F(/0) _ ae -W))) sin(3(F(/i)) 

= -fc S(fB(ft)) siii(3(F B (li))) 1 (8) 

where F(h) = ip(t = hM). General solution to the above equation is complicated if not 
impossible to obtain. Therefore we will consider a special case which is useful in analyzing 
the ensuing example. We will consider the case where a and c phases are symmetric and b 
phase is also symmetric with respect to the midpoint between the two. Furthermore we will 
consider the scaling limit where corrections in 0(L~ 2d ) is neglected. In that case we have 

cosh(3?(F(/i))) cos(9?(F(/i))) = -6/2, (9) 
sinh(3?(F(/i))) sin(3(F(/i))) = 0. (10) 

If b/2 < 1, we have zeros as = cos" 1 (-6/2)/L a! /i and U(h) = while for 6/2 > 1, 

^s(h) = (7r(l + 2k))/L d [i and 9ft(/i) = ± | cosh -1 6/2 | /L d jjL. This means that crossing the 
triple point from two phase region(6/2 < 1) to the single phase region(6/2 > 1), zeros on 
the unit circle disappear. Only for small sized systems and close to the triple point, we see 
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two circles separating out from the unit circle which are sort of ghost circles that separate 
three phases, two unstable symmetric and the new stable phases. 

We have tested the above conclusion explicitly by calculating zeros of the partition 
function of the Blume-Capel model. Consider the energy of the Blume-Capel model [7||§, 
E = —H + D J2i Sf — J J2<i,j> SiSj of N spins on a square lattice of side L(N = L 2 ) 
on a torus. H, D and J are usual external parameters and Si are spin variables which take 
three values, ±1 and 0. If we designate by N + and iV_ the number of +1 spins and — 1 
spins, and N e = [{q/2) J2i Sf — J2<i,j> SiSj], we can write the energy E = —H(N + — iV_) + 
A((iV + + AL) + JN e where q is the coordination number and A = D — J q/2. We have 
defined energy N e such a way to make the ground state a triple point, in the absence of 
external fields H and A. The partition function can be calculated according to 

Z{T,H,A)= ]T e- E ^ N +' N -^ kBT cu(N e ,N + ,N.), (11) 

N e ,N+,N- 

where w(N e , N + , N_) is the number of configurations of given N e ,N + and iV_. We have 
calculated exact u(N e , N + , A r _)'s for all N + , N_ for N e = 0, 1, • • • , 6. These cj's are sufficient 
to calculate the partition functions at temperatures up to '/ J = 0.2 for sizes up to L = 70 
with the machine precision. The phase diagram in T — H — A space looks like Fig. [T| [[J . 

Defining order parameters, iV = iV + + iV_ and M = N + — AL and their conjugate 
field variables, ( = j3(A — A G ) and r) = (3{H — H a ), we can write the partition function 
as Z((,rj) = J^n m e~ < ' N+vM Z (N, M). In the above H Q and A Q are the values used to 
calculate the microcanonical partition function Z a (N, M) which is defined as Z D (N, M) = 
e -0A o N+/3H o M ^ Ne e -/ 3J7V ^(iV e , (N + M)/2, (N - M) j2). If H and A are chosen to be the 
transition point, ( and r\ will vanish at the transition point. 

In order to maximize the finite size effect, we demonstrate zeros at the highest allowable 
temperature, (3 = 5.0 where the high-order cumulants are appreciable while the contributions 
to Z a (N,M) from terms of N e higher than 6 are still negligible. The transition points are 
determined by equating the free energies of the coexisting phases for the system size L ~ 60. 
We found the triple point at this temperature, Pt = (Co = /^A = —0.0000454657, r\ = 
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(3H = 0)(Fig. 1(e)). A typical asymmetric transition point which we chose for the demon- 
stration is P A = (Co = l,Vo = 1.0000061392784119), a point like Fig. 1(f). 

In Fig. 2. we display ln(Z G (iV + + AL, N + - JV_) by the size of the filled circle in (N+, AL) 
plane at various points in the phase diagram. It should be noted that in these figures M 
and N vary along the two diagonal lines. 

We found that zeros in the complex z v (= e' 7 )-plane at phase points like Fig. 1(a) along 
the r] = line up to the triple point are distributed more or less uniformly on the unit 
circle for all L as predicted by the theorem. On the other hand there are no zeros along 
the zero-field line in z^(= e^-plane in the direction parallel to rj = line as long as the 
phase points remain far away from the triple point. At Pt we found zeros in 2^-plane in the 
direction indicated by Fig. 1(e). form a circle of varying radius which is dependent on the 
size of the system as shown in Fig. 3A. Since the transition at Pt in the direction parallel 
to r] = line is asymmetric, it might be concluded that the partition functions of the two 
phases are added up with unequal prefactors. However the asymmetry in this case is due 
to the fact that it is a transition from the two-phase region to the single phase region, and 
the asymmetry factor a in Eq.(^) is indeed exactly 2 if we fit the concentric circles by the 
predicted relations, r = exp[ln(a)/(2L d /i)]. 

On the other hand in the z^-plane the uniformly distributed 2N zeros pair up and 
converge to each other as ( Q approach the triple point and they eventually double up before 
they bifurcate into two sets of concentric circles of N zeros each at Pp. Fig. 3B is the plot 
of zeros at P T = (£ = 0.01088, r\ = 0), a point slightly to the right of Fig. 1(b). The size 
dependence in this plot is not from the L's in F(h) in Eqs.(|^) and ([II]) which would give 
converging circles rather than divergence. The origin of the divergence is the L-dependence 
of the the parameter b of Eq.([|). Since the phase point P' T is far enough from the triple point 
so that b(L) oc exp((( (P T ) - ( (P T ))L d ), r diverges as r oc exp((( (P T ) - Co(ir))^ d )- Since 
circles of both sets converge to the unit circles if we move r] from r\ = 0, to a coexistence 
point on either B+ or B— plane of Fig. 1, the two sets are actually finite-size ghost circles 
that separate three phases bounded by B+ and B— plane. 
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Fig. 3C is a plot of zeros in z^-plane at Pa and at two nearby points, P^ = (Co = 
Co (-Pa) ± 0.000005,^ o = t) (Pa)) Although zeros for only two sizes, L = 8 and L = 40 are 
shown, zeros of all sizes collapse into single curves. This proves that a — 1, which implies 
that the partition functions of two phases are added up without a prefactor. 

Fig. 3D is the plot of exact zeros at Pa, together with approximate zeros calculated using 
Eqs.(f|) and (§]). We calculated ji for I = 1 through 6 and approximate zeros of successive 
order which include 7/ of successively higher orders. The reason that we see only 4 sets of 
approximate zeros is because a = 1 or R = makes the third and fifth terms of Eq.(f|) 
vanish. It should be noted that R = makes even-order terms of Eq.(|5|) vanish instead. 
This again proves a = 1. It also shows that the locus of zeros for large 9 is affected more and 
more by the discontinuities of higher order derivatives of the free energy, 7V s as predicted 
by the theorem. 

In conclusion we have shown that the extended circle theorem for first-order transitions 
follows from a very simple property that the distribution function for the order parameter 
is multiply peaked at the coexistence point. In fact zeros of the partition function can be 
expressed in terms of the position and shape of these peaks through cumulants of each phase. 
We have also found the finite-size effect on the zeros of the partition function which can be 
useful in determining the order of transition as well as the accurate estimate of the transition 
point from the numerical data of finite-size systems which can be easily obtainable from a 
Monte Carlo technique of Ref. |nj for example. Finally we have found that the partition 
functions are added up without a prefactor in the case of asymmetric phase transition in 
the Blume-Capel model, at least at low temperatures. 
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FIGURES 

FIG. 1. Schematic phase diagram of Blume-Capel model. The symmetric coexistence plane 
A lies in T — A plane bounded by the triple line OP and the critical line PR. Two asymmetric 
coexistence planes B+ and B— bifurcate from the triple line symmetrically. P is the tricritical 
point. 

FIG. 2. \n(Z a {(N + + N_,N+ - JV_)) on (N+,N-) plane for L = 16. I. At Co = -l,ffo = 0, a 
symmetric two-phase transition point corresponding to Fig. 1(a). II. At Pr- III. At Co = F?7o = 0) 
a symmetric single phase point far away from the triple point corresponding to Fig. 1(c). IV. At 

Pa- 
fig. 3. Zeros of partition function. Circles of a solid line are unit circles in all 4 figures. A. 
6 sets of N(= L 2 ) zeros for the system size, L = 8, 10, 12, 14, 16 and L = 60 in decreasing size of 
dots at Pt in z^-plane. The radial distance from the unit circle is magnified 50 times. Even at 
this magnification zeros for L = 60 is too close to the unit circle to be discernible. B. Zeros P' T 
in z^-plane for L = 8, 10, 12 and L = 14. Distance from the unit circle for outer circles is scaled 
down by a factor of 1/25 to bring them into the plotting area. C. Zeros in z^-plane at Pa (the circle 
tangent to the unit circle), P^"(the one inside of the unit circle) and (the largest circle which 
crosses the unit circle). In C and D the radial distance from the unit circle are again magnified 
this time 50000 times! D. Exact and approximate zeros at Pa- Zeros for two sizes, L = 8 and 40 
are shown. The unit circle is the zeroth order approximation, the innermost is the second order 
approximation, zeros closest to the unit circle is the fourth, and the sixth is closest to the exact 
zeros, the circle in the middle. Similar plots for zeros in z^-plane (in the direction of Fig. 1(d)) are 
obtained at this point. 
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